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ABSTRACT 

We discuss, using simple analytical models and MHD simulations, the origin and pa- 
rameters of turbulence and magnetic fields in galaxy clusters. Any pre-existing tangled 
magnetic field must decay in a few hundred million years by generating gas motions 
even if the electric conductivity of the intracluster gas is high. We argue that tur- 
bulent motions can be maintained in the intracluster gas and its dynamo action can 
prevent such a decay and amplify a random seed magnetic field by a net factor typi- 
cally 10 in 5 Gyr. Three physically distinct regimes can be identified in the evolution 
of turbulence and magnetic field in galaxy clusters. Firstly, the fluctuation dynamo 
will produce microgauss-strong, random magnetic fields during the epoch of cluster 
formation and major mergers. At this stage pervasive turbulent flows with r.m.s. ve- 
locity of about 300 km s -1 can be maintained at scales 100-200 kpc. The magnetic 
field is intermittent, has a smaller scale of 20-30 kpc and average strength of 2/iG. 
Secondly, turbulence will decay after the end of the major merger epoch; we discuss 
the dynamics of the decaying turbulence and the behavior of magnetic field in it. 
Magnetic field and turbulent speed undergo a power-law decay, decreasing by a fac- 
tor of two during this stage, whereas their scales increase by about the same factor. 
Thirdly, smaller-mass subclusters and cluster galaxies will produce turbulent wakes 
where magnetic fields will be generated as well. Although the wakes plausibly occupy 
only a small fraction of the cluster volume, we show that their area covering factor 
can be close to unity, and thus they can produce some of the signatures of turbulence 
along virtually all lines of sight. The latter could potentially allow one to reconcile 
the possibility of turbulence with ordered filamentary gas structures, as in the Perseus 
cluster. The turbulent speeds and magnetic fields in the wakes are estimated to be 
of order 300 km s -1 and 2/xG, respectively, whereas the turbulent scales are of order 
200 kpc for wakes behind subclusters of a mass 3 x 10 13 Af Q and about 10 kpc in the 
galactic wakes. Magnetic field in the wakes is intermittent and has the scale of about 
30 kpc and 1 kpc in the subcluster and galactic wakes, respectively. Random Fara- 
day rotation measure is estimated to be typically 100-200 radm~ 2 , in agreement with 
observations. We predict detectable polarization of synchrotron emission from cluster 
radio halos at wavelengths 3-6 cm, if observed at sufficiently high resolution. 
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1 INTRODUCTION 

Intergalactic gas in clusters of galaxies is magnetized (see 
reviews by Kronberg 1994; Carilli & Taylor 2002; Govoni 
& Feretti 2004). The number of clusters that exhibit de- 
tectable synchrotron emission is relatively small, but it is 
believed that a magnetic field is present in most clusters 
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(unlike relativistic electrons). Therefore, a more informa- 
tive observational tracer of intracluster magnetic fields is 
Faraday rotation of the polarization plane of background 
radio sources and central radio galaxies. Clarke, Kronberg 
& Bohringer (2001) conclude, from their statistical study of 
Faraday rotation in 16 galaxy clusters, that random mag- 
netic fields of a strength (5-10) x (Z/10 kpc) ~ 1/2 fiG (with I 
the field scale) permeate the gas within about 500 kpc from 
the cluster centre; the area covering factor of magnetic fields 
is close to unity. The Faraday rotation measures detected 
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are about 200radm 2 for the lines of sight through clus- 
ters' central parts and about lOOradm -2 farther out (see 
also Clarke 2004; Johnston-Hollitt & Ekers 2004). Important 
constraints on intracluster magnetic fields come from lim- 
its on X-ray emission produced from microwave background 
photons by inverse Compton scattering off relativistic elec- 
trons (Bagchi, Pislar & Lima Neto 1998; Sarazin 1988). Ob- 
servational evidence is compatible with a random magnetic 
field of r.m.s. strength of 1-10 fj,G and coherence length of 
about 10-20 kpc. 

Faraday rotation maps of a number of radio galaxies 
in clusters have also been analyzed. Eilek and Owen (2002) 
studied Faraday rotation maps of radio sources in the centres 
of the Abell clusters A400 and A2634, and found patches of 
RM fluctuations on scales about 10-20 kpc. Assuming this 
to be also the coherence scale for the field, these authors 
deduce field strengths of 1-4 pG. Vogt and Enfilin (2003, 
2005), using a novel technique (Enfilin & Vogt 2003), es- 
timate magnetic field strength to be 3 [iG in A2634, 6 pG 
in A400 and 7 fiG in Hydra A. They obtain field correla- 
tion length of 4.9 kpc, 3.6 kpc and 3 kpc for these clusters, 
respectively. It cannot be excluded, however, that the Fara- 
day rotation of cluster radio sources is contaminated by that 
arising in dense turbulent cocoons around the radio galax- 
ies, rather than in the intracluster medium proper (Rudnick 
and Blundell 2003; see Enfilin et al. 2003 for another view). 
For this reason, the statistical studies of Faraday rotation 
of background radio sources, referred to above, provide per- 
haps a more convincing evidence for cluster-wide magnetic 
fields and their properties. 

The origin of the cluster magnetic fields remains un- 
clear. Carilli & Taylor (2002) (see also Tribble 1993a) argue 
that the small value of electric resistivity of the intraclus- 
ter plasma guarantees that the decay time of magnetic field 
will be comparable to or exceed the cluster lifetime. They 
conclude that any magnetic field (e.g., that captured by a 
cluster during its formation) would survive for a long time. 
However, any inhomogeneous magnetic field will drive mo- 
tions via the Lorentz force, and the motions will decay, plau- 
sibly in the form of decaying MHD turbulence (e.g., Biskamp 
2003; see below). The turbulent decay time is comparable to 
the eddy turnover time of the largest eddies, about 10 8 yr, 
irrespective of the resistivity or viscosity of the gas. (Here 
we have taken a coherence scale of order 10 kpc for the field 
and the induced motions, with associated turbulent veloc- 
ities of order 100 km s -1 ). Although the energy density in 
MHD turbulence decays with time as a power law, this time 
scale is still much shorter than the typical age of a cluster, 
which is thought to be several billion years. Therefore, one 
has to provide explicit explanation of the origin and persis- 
tence of magnetic fields in the clusters; reference to the low 
Ohmic resistivity of the intracluster plasma is not sufficient 
if the gas is turbulent or the magnetic field is tangled. 

An obvious option to explain intergalactic magnetic 
fields is to consider magnetic fields stripped from galaxies. 
Since the intracluster gas is enriched with metals, at least 
part of it originates in galaxies (Sarazin 1988). However, the 
strength of magnetic field produced by the stripping cannot 
exceed ~ 0.1 /iG even in the cores of rich clusters and even 
if spiral galaxies with relatively strong large-scale field are 
involved (see Appendix 0. Another possibility is that the 
intracluster field is supplied by active galaxies within the 



cluster. As we discuss in Appendix^ this mechanisms can 
provide relatively strong magnetic field but fails to explain 
how the field can be maintained against turbulent decay. 
In addition, it is not quite clear how efficiently the mag- 
netized relativistic plasma of the radio lobes can be mixed 
with the thermal intergalactic plasma and what would be 
the resulting scale of magnetic field. Altogether, the above 
mechanisms can only provide suitable seed magnetic field 
for the dynamo action in the intracluster plasma. 

In most astrophysical systems, like disc galaxies, stars 
and planets, rotation is crucial for maintaining their mag- 
netic fields, both by providing strong shear and by making 
(when coupled with stratification) random flows helical, and 
hence leading to mean-field dynamo action. However, galaxy 
clusters are believed to have fairly weak rotation (if any at 
all), so one has to appeal to some other mechanism for un- 
derstanding cluster magnetism. 

Another possibility to generate magnetic fields is related 
to the fluctuation dynamo action (Batchelor 1950; Kazant- 
sev 1967) , where random flow of electrically conducting fluid 
generates random magnetic field (Zeldovich, Ruzmaikin & 
Sokoloff 1990). This mechanism does not require any rota- 
tion or density stratification, and only relies on the random 
nature of the flow; so fluctuation dynamos can be active 
in virtually any turbulent environment where the plasma 
is ionized. Apart from the randomness of the flow, it is re- 
quired that the magnetic Reynolds number is large enough, 
i.e., that the electric conductivity is high enough, and/or 
plasma motions are sufficiently intense, and/or their scale is 
sufficiently large. 

The earliest theories of intracluster magnetic fields were 
based, implicitly or explicitly, on fluctuation dynamo theory 
under the assumption that galaxy clusters are steady-state 
turbulent systems (Jaffe 1980; Roland 1981; Ruzmaikin, 
Sokoloff & Shukurov 1989; Goldman & Rephaeli 1991; De 
Young 1992). The source of turbulence adopted by several 
authors were turbulent wakes of the cluster galaxies. This 
picture has to be reconsidered for several reasons. 

Firstly, turbulence from galactic wakes can fill the clus- 
ter volume only if the effective galactic radius is of order 
10 kpc (e.g., Ruzmaikin et al. 1989), i.e., if the interstellar 
gas is not stripped by the ram pressure of the intraclus- 
ter gas. If the gas stripping is complete, the wake is only 
produced by gravitational accretion (Bondi 1952), and its 
radius is about the accretion radius r g = 2GM/(c 2 + V 2 ) 
where c s is the speed of sound, V is the galactic speed, M 
is the galactic mass, and G is Newton's gravitational con- 
stant. For c s = 10 3 kms~ 1 , V « Cs and M = lO n M , the 
gravitational accretion radius r g ~ 0.5 kpc is much smaller 
than both the galactic radius and the apparent scale of the 
random magnetic field in the intergalactic gas (both usually 
assumed to be of order 10 kpc). Hence, sufficiently strong, 
volume-filling turbulent wakes whose width is comparable 
to the galactic size can only arise if the galaxies retain sig- 
nificant amounts of their interstellar gas. If the stripping of 
interstellar gas by ram pressure is efficient, galactic wakes 
are rather weak (Portnoy, Pistinner & Shaviv 1993; Balsara, 
Livio & O'Dea 1994; Acreman et al. 2003 and references 
therein; see however Toniazzo & Schindler 2001 who argue 
that the stripping efficiency is exaggerated in the above pa- 
pers). Therefore turbulence generated in galactic wakes may 
not fill the cluster volume (see also Sect. l2.3~H . 
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Further, numerical simulations of De Young (1992), us- 
ing a closure model, gave pessimistic estimates for mag- 
netic fields produced by the dynamo when the turbulence 
is induced by galactic wakes. However, it is not clear if the 
resolution of those simulations (i.e., the effective magnetic 
Reynolds number) was high enough to obtain dynamo ac- 
tion, so this objection to the dynamo models is question- 
able. Recent direct simulations of dynamo action in tur- 
bulent flows (Haugen, Brandenburg & Dobler 2003, 2004; 
Schekochihin et al. 2004) have confirmed the efficiency of 
dynamo action in random non-helical flows. 

Perhaps more importantly, cluster dynamics has been 
reconsidered recently, and the emerging picture is very dif- 
ferent in that clusters may not be relaxed systems, but still 
remain in the state of formation via major mergers and ac- 
cretion of smaller-mass subclusters. Numerical simulations 
and recent observations strongly suggest that a random flow, 
perhaps of turbulent nature, can be maintained for a few 
crossing times of the forming galaxy cluster (Norman & 
Bryan 1999). Roettiger et al. (1999a,b) have found that mag- 
netic field can be amplified by these motions (see also re- 
sults obtained with SPH simulations by Dolag, Bartelmann 
& Lesch 1999, 2002). However, the resolution of the simu- 
lations is still poor, and quantitative estimates of the tur- 
bulence parameters and especially of its effects on magnetic 
field are very uncertain. 

In addition to the volume-filling flow produced during 
the cluster formation, significant random flows can still be 
generated in wakes behind infalling subclusters and cluster 
galaxies. These are not expected to fill the cluster volume, 
but we argue below that they can have significant area cov- 
ering factor. There is also a possibility that radio galaxies 
can stir the intracluster gas as their plasma buoyantly rises 
through the gas (Briiggen et al. 2002; Enfilin & Heinz 2002). 
Another possible consequence of radio galaxy jets and/or 
lobes propagating at subrelativistic speeds through the clus- 
ter plasma, is the generation of turbulence in a cocoon sur- 
rounding the radio source (see for example Reynolds, Heinz 
& Begelman 2002). 

The content of the paper is as follows. We consider the 
evolution of turbulence in the intracluster gas of a galaxy 
cluster during and after its formation in Sect. |5] Random 
flows produced during the merger epoch can lead to mag- 
netic field generation via the fluctuation dynamo as dis- 
cussed in Sect. In Sects l2~2l and HI we present evidence 
that the decay of both magnetic and kinetic energies after 
the epoch of major mergers will be a power law in time, 
rather than exponential, because of the turbulent nature of 
the flow. During the decay phase of the turbulence, the cor- 
relation scale of the magnetic field will grow (Frisch 1995; 
Olesen 1997; Biskamp & Miiller 1999; Christensson, Hind- 
marsh & Brandenburg 2001), which slows down the decay of 
the Faraday rotation measure produced in the intracluster 
gas (Sect. l5~Tll . In Sect. 12.31 we argue that significant turbu- 
lence and magnetic field amplification can occur in turbulent 
wakes of smaller-mass subclusters and cluster galaxies. The 
turbulent wakes may not fill the volume, but can cover the 
cluster's projected area (Sect. |2~3.21 . We present order-of- 
magnitude estimates of the parameters of the random flows 
at various stages of the cluster evolution in Sect. |5] and 
of magnetic fields generated by the flows, in Sect. [31 which 
are further substantiated by numerical simulations of dy- 



namo action in driven and decaying random flows discussed 
in Sect. [I] Faraday rotation measure and polarized radio 
emission produced by these magnetic fields are discussed in 
Sect. |S| Our results are summarized in Sect. [U and Tabled 



2 INTRACLUSTER TURBULENCE 

An upper limit on the turbulent velocity in a steady-state 
galaxy cluster follows from the requirement that the rate of 
dissipation of the turbulent energy should not exceed the 
X-ray luminosity of the cluster Lx, i.e., ^Vq/Io < Lx/M g , 
where vo and Zo are the turbulent speed and scale, respec- 
tively, and M g is the mass of the intergalactic gas. This 
yields 

Wo <l 80 ^(_i2_j ( 104 f^ g/s ) (yq^) • W 

This restriction applies to a steady state, and stronger tur- 
bulence can be driven in an evolving cluster, where the en- 
ergy released by the decay of turbulent motions heats up 
the gas. Therefore, turbulent velocities significantly exceed- 
ing the above value can be considered as an indication of 
the cluster's ongoing evolution. The turbulent nature of the 
flow is important, however: Eq. Q only applies to turbulent 
flows and has to be reconsidered in the case of a random flow 
without turbulent energy cascade. As we argue below, tur- 
bulent motions during the epoch of cluster formation are in- 
tense enough to violate the constraint Q because they keep 
evolving and their dissipation contributes to the heating of 
the intracluster gas to the virial temperature. 

We now discuss various sources of turbulence in galaxy 
clusters. 

2.1 Turbulence produced during cluster formation 

Theories of hierarchical structure formation suggest that 
clusters of galaxies have been assembled relatively recently. 
iV-body simulations indicate that the clusters form at the 
intersection of dark matter filaments in the large-scale struc- 
ture, and result from both major mergers of objects of com- 
parable mass (of order 1O 15 M0) and the accretion of smaller 
clumps onto massive protoclusters. It is likely that intense 
random vortical flows, if not turbulence, are produced in the 
merger events (Kulsrud et al. 1997; Norman & Bryan 1999; 
Ricker & Sarazin 2001). These would originate not only due 
to vorticity generation in oblique accretion shocks and insta- 
bilities during the cluster formation, but also in the wakes 
of the smaller clumps. In this section we summarize some of 
the work on the cluster-wide turbulence resulting from ma- 
jor mergers. Its consequences for the generation of cluster 
magnetic fields are considered in Sect. [3] 

There have been several numerical simulations of gas 
dynamics during the formation of galaxy clusters. Norman 
& Bryan (1999) find that the intracluster medium becomes 
turbulent during cluster formation, with turbulent velocities 
of about 400 km s _1 within 1 Mpc from the centre of a cluster 
and eddy sizes ranging from 50 to 500 kpc; the random flow 
is volume filling (see also Sunyaev, Norman & Bryan 2003). 
In the cluster merger model of Ricker & Sarazin (2001), ram 
pressure displaces gas in the cluster core from the bottom of 
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the potential well. The resulting convective plumes produce 
large-scale disordered motions with eddy size up to several 
hundred kiloparsecs; even after 15 Gyr of evolution, turbu- 
lent velocities in the inner parts remain at a surprisingly high 
level of 10-40% of the sound speed (i.e., 100-400 kms" 1 ), 
and grow up to the sound speed in the outer parts. 

Observational evidence of intracluster turbulence is 
scarce. From analysis of pressure fluctuations as revealed in 
X-ray observations, Schueker et al. (2004) argue that the in- 
tegral turbulent scale in the Coma cluster is close to 100 kpc, 
and they assume a turbulent speed of 250 km s -1 at that 
scale. The nonthermal broadening, by the turbulence, of X- 
ray spectral lines of ionized iron can be detectable with fu- 
ture X-ray observatories (Inogamov & Sunyaev 2003). 

The spatial resolution of all available simulations is 
rather coarse (of order 10 kpc or worse) and only compa- 
rable with the apparent scale of the intracluster magnetic 
field. Even if the random nature of the resulting flow is ob- 
vious, it is not clear if it will evolve into developed turbu- 
lence, i.e., if the turbulent cascade is established to carry 
energy to small scales where it dissipates into heat. From 
the viewpoint of magnetic field generation, the presence of 
turbulence as such is not required; the randomness of the 
flow is sufficient (Kazantsev 1967; Zeldovich et al. 1990). 
However, the dynamics of the flow and its evolution do de- 
pend on whether or not the turbulent cascade persists (see 
Sect. . 

The flow can become turbulent if the kinematic 
Reynolds number in the intracluster gas, Re, is large enough. 
Following Sarazin (1988) and Ricker & Sarazin (2001), an 
estimate of Re can be obtained as 

Re= ~3 = 3X 2 

V C S A<3 AO 

where v = gC a Ao~ is the effective kinematic viscosity, c s is the 
speed of sound (assumed to be close to the thermal speed), 
M is the Mach number, A is the ion mean free path, and 
subscript '0' refers to the energy-range values. Here we have 
introduced parameter S that quantifies the poorly under- 
stood behaviour of viscosity in the intracluster plasma. The 
standard, Spitzer's value of A can be written as 

A - 5kpC (lO3kms-0 (lO-^cm-O ' 
with n the electron number density. For J\A ~ 1, lo ~ 100- 
500 kpc and A> lkpc, this yields Re < (300-1500)5- 1 . This 
estimate is, however, suspect because the mean free path 
is comparable to the scale of inhomogeneities in the gas. 
It is also clear that even a weak seed intracluster magnetic 
field could strongly reduce the effective viscosity and make 
it anisotropic. The effective Reynolds number can be sig- 
nificantly larger if any shorter length scale plays the role 
rather than the Coulomb mean free path, or a frequency 
higher than the ion collision frequency (these may be asso- 
ciated with plasma instabilities and/or waves). Schekochihin 
et al. (2005b) argue that the firehose and/or mirror insta- 
bilities can provide the effective diffusion in the magnetized 
plasma; the corresponding length scale is the ion gyroradius, 
a quantity normally much smaller than the mean free path. 
These uncertainties are allowed for, in a heuristic manner, 
by choosing 8 < 1; Fabian et al. (2005) scale their results by 
S — 0.1. This prescription can be an oversimplification as it 
omits plausibly important physical effects arising from the 
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Figure 1. The evolution of the spectrum of decaying turbulence, 
where the spectral exponent remain constant both in the inertial 
range, where it is equal to —5/3, and at large scales where it 
is equal to s. Solid, dashed and dotted lines show the spectrum 
at consecutive times. As the energy dissipates, the energy-range 
wave number kg reduces. 

anisotropy of viscosity in magnetized plasma (cf. Schekochi- 
hin et al. 2005b). However, even if one cannot readily provide 
a confident estimate of Re, Eq. suggests that it will be 
large enough as to ensure that random motions driven by 
major merger events can become turbulent. Only this fact is 
important for our purposes in this paper. Our simulations of 
a flow driven by external random force (Sect. have been 
performed for Re = 100-400, and they do indeed show a flow 
with broad range of scales typical of turbulent flows (and so 
inertia forces dominate over viscosity), even if the inertial 
range is not wide at these modest Reynolds numbers. 

To summarize, the above results seem to converge to the 
following picture. Random motions driven in major merger 
events have the typical initial speed of i>oi — 300kms _1 
and scale Zoi ~ 100-200 kpc, so that the turnover time of 
the energy-range eddy is toi = 0.3-0.6 Gyr. The random 
motions will be maintained at this level during the major 
merger epoch, whose duration can be as large as tf ~ 3- 
5 Gyr. (The notation is motivated in Sect. I2.2T ) 

2.2 Dynamics of decaying turbulence 

The random flows produced by major mergers will not re- 
main statistically steady after the end of the merger event. 
Unlike a laminar flow that decays exponentially in time due 
to viscosity, turbulent kinetic energy decays slower, as a 
power law (e.g., Landau & Lifshitz 1975, Frisch 1995). The 
reason for this is that kinetic energy mainly decays at small 
scales, to where it is constantly supplied by the turbulent 
cascade. As a result, the energy decay rate depends nonlin- 
early on the energy itself, which makes the decay a power law 
in time. (The corresponding calculation is provided below.) 
Our simulations of Sect.[I]confirm that the power-law decay 
occurs even for the Reynolds number as small as Re ~ 100. 

In this section we review simple models of decaying hy- 
drodynamic turbulence. The effects of magnetic field on the 
flow can be neglected at early stages when magnetic field is 
still weak. In Sect.^Jwe present numerical simulations where 
the effects of magnetic field on the flow are fully allowed for. 

Consider an initial spectrum of turbulence shown with 
solid line in Fig.Q where Ej, = Ck s (with s > — 1) at scales 
k < fco, and Ek oc k~ 5 ^ 3 at smaller scales as in Kolmogorov 
turbulence. Here kEu is the specific energy per unit loga- 
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rithmic interval in the fc-space. It is related to the turbulent 
velocity at wave number k via Vk oc (kE k ) l/2 ; E k has a max- 
imum at a certain wave number ho, which is therefore called 
the energy-range wave number. Turbulent flow remains in a 
statistically steady state despite viscous dissipation at small 
scales if it is driven at larger scales. When such a driving 
ceases, the turbulence decays. 

Motions at small scales are the first to be affected by 
viscosity. It is reasonable to expect that the exponent s and 
constant C are preserved during the decay (e.g., Sect. 7.7 
of Frisch 1995). Consider the time-dependent total spe- 
cific turbulent energy E(t), which is approximately equal 
to that at the energy-range scale, E(t) ~ ^v 2 , if the inertial 
part of the spectrum is steep enough. On the other hand, 
E(t) = E k dk oc fcg +1 . As long as the Reynolds number 
remains large, and so viscosity at ko is negligible, energy at 
the energy-range wave number ko(t) mainly decays because 
it cascades to smaller scales. Hence, the energy loss rate is 
given by dE/dt = —Vg/to ~ — vfao, where to m (voh))^ 1 . 
Since v oc E 1/2 and k oc _E 1 /( s+1 ) j the evolution of the 
total turbulent energy is governed by 
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(3s+5)/[2(s+l)] 
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where A is a certain constant related to C. This equation 
can easily be integrated (for constant C) to yield (see Fig.0 
asymptotic forms (applicable at large t) for the energy decay 
law 

E{t) oc r a , 

and for the energy-range scale growth law 

lo oc kg 1 oc t 13 , 

where 

,s + l „ 2 



a = 2- 



s + 3 



s + 3 



The value s — 2 gives a 'white noise' spectrum at large 
scales, where the three-dimensional spectrum k~ 2 E k is fiat. 
In this case 



a = 6/5 



= 2/5 



Another possibility often considered in this context is s = 4, 
where one gets a = 10/7 and (3 = 2/7 (Skrbek & Stalp 2000; 
Touil, Bertoglio & Shao 2002). 

There are arguments suggesting that s — 2 is a better 
acceptable value than s = 4; in particular, the coefficient C 
is time-dependent for s — 4, and this makes the energy decay 
significantly slower than that obtained above (Frisch 1995) . 
Also s = 4 would be relevant for an incompressible flow, 
whereas the cluster flows are expected to be mildly com- 
pressible, and s — 2 could be more appropriate. The turbu- 
lence decay is sensitive to the detailed physical nature of the 
system, and it is often slower than derived above. The decay 
of the MHD turbulence can be significantly slowed down if 
the system has non-zero invariants such as magnetic helic- 
ity and/or cross-helicity (Biskamp 2003). If the intracluster 
seed magnetic fields are due to stripping of the galactic mag- 
netic fields, then they may have both types of helicities. The 
decay law is also sensitive to the relation between the turbu- 
lent energy-range scale and the size of the system; the decay 
speeds up to become E oc t~ 2 when the two scales become 



comparable and the value of ko cannot decrease any further 
(Skrbek & Stalp 2000; Touil et al. 2002). 

We adopt a — 6/5 and j3 = 2/5 for numerical estimates 
in what follows, so that 
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for t — t{ 3> £oi , where subscript T refers to the start of the 
evolution, toi is a certain dynamical time scale, which can 
be identified with the initial turnover time of the energy- 
containing eddies, toi = Zoi/«oi, subscript '0' refers to the 
energy-range (correlation) scale of the motion (which varies 
with time), and the decay starts at time t = tt when the 
flow forcing ceases. This decay is faster than in many other 
models of decaying turbulence; thus, our conclusions will be 
rather conservative with respect to the intensity of turbu- 
lence at late times. 

With the above decay law of turbulence, the Reynolds 
number evolves slowly as 



R e oc 



-1/5 



t - tf > t I 



Allowing for the initial period tf = 3 Gyr of sustained tur- 
bulence, Re decreases only by a factor of 1.4 after the to- 
tal evolution time of t = 6 Gyr for loi = 150 kpc and 
Voi = 300 kms -1 (yielding ioi = 0.5 Gyr). 



2.3 Minor mergers and turbulent wakes 

Consider the infall of relatively small subclusters of mass 
m into an already formed cluster of mass M. Define 
d 2 p/(d(lnm) dt) as the probability that in a time dt a sub- 
cluster, whose mass belongs to a logarithmic mass interval 
[mm, lnm + d(lnm)], merges with the bigger cluster of mass 
M. The merger rate d 2 p/(d(ln m)dt) scales with the subclus- 
ter mass as (Lacey & Cole 1993) 



d p _i/2 
oc m ' 



d(ln m) dt 



for — < 1 
M 



(•») 



Thus, the merger rate of masses of order 1O 13 M0 is about 
10 times larger than that for 10 15 M©. If major mergers of 
masses of order 10 15 A/q occur once in 3 Gyr, the time in- 
terval between mergers with 1O 13 M0 subclusters will then 
be of the order of 0.3 Gyr (see also Norman & Bryan 1999). 
Such minor mergers are thought to play an important role 
in explaining the observed cold fronts in clusters (Heinz et 
al. 2003; Motl et al. 2004 and references therein). They can 
also generate turbulence in the wake of a moving subcluster 
(cf. recent simulations of Takizawa 2005). Turbulence gen- 
erated by subclusters was suggested by Norman and Bryan 
(1999) to be a major source of the random motions observed 
in their simulations of cluster formation. Here we examine 
this issue further with analytical estimates, using parame- 
ters of galaxy clusters and smaller structures obtained from 
hierarchical theories of structure formation. 



2. 3. 1 Ram pressure stripping 

The subclusters contain gas which can be partially stripped 
by hydrodynamic interaction with the cluster gas (by 
ram pressure stripping and via hydrodynamic instabilities) 
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(Fabian & Daines 1991; Acreman et al. 2003). A simple cri- 
terion for the radius Ro within which the subcluster gas 
remains unstripped can be obtained as follows. The ram 
pressure force exerted on a gas sphere of radius Ro is equal 
to p c v 2 c ivRo, where p c is the intracluster gas density and v BC 
is the speed at which the subcluster moves through the in- 
tracluster gas. Following Fabian & Daines (1991), we note 
that the gravitational restoring force per unit area due to 
the subcluster mass is comparable to the gas pressure in the 
subcluster, assuming that the subcluster is in hydrostatic 
equilibrium. The gas sphere will be removed from the sub- 
cluster if the restoring force is smaller than the ram pressure 
force. Thus, a local criterion for retaining the gas at a dis- 
tance Ro from the subcluster centre is 



PcV sc < fp sc (R )u 



(6) 



where p S c(R) is the gas density distribution of the subcluster, 
R is the subcluster's spherical radius, / is a numerical factor 
of order unity, and u is the gas velocity dispersion within the 
subcluster. We adopt, for illustrative purposes, gas density 
profiles for the cluster and subcluster, respectively, of the 
form 



Pc 



[1 + (r/r c ) 2 ] 



p a c(R) = 



[1 + {R/R sc ) 2 ] 



where p CQ and p sco are the respective central gas densities 
and r c and R sc are the corresponding gas core radii. (These 
correspond to the standard /3-profile with the slope param- 
eter P = 2/3 - Sarazin 1988). From Eq. ®, the subcluster 
gas is retained at radii smaller that Ro, where 



' Ro\ 2 = f PscpV 2 

\ Rsc J PcQ^SC 



- 1 



(T) 



Takeda, Nulsen & Fabian (1984) suggest / ~ 2. 

Parameters of clusters and subclusters that enter 
Eq. Q vary broadly in both observed and simulated clus- 
ters. Suitable values can be selected as follows. For example, 
consider subclusters predicted by the hierarchical theory of 
structure formation (Peebles 1980; Padmanabhan & Subra- 
manian 1992; Padmanabhan 1993). Suppose that the initial 
density fluctuations can be described as a Gaussian ran- 
dom field with the r.m.s. density contrast a m (m), where 
m is the mass of the structure. In the hierarchical theory, 
(r m (m) oc m _ ' 3+n ' / ' 3 with n close to —1 at the cluster scales 
and to —2 at the galactic scales. For a density fluctuation 
which is p times the above r.m.s. value, the following scal- 
ing laws can be obtained: r v j r oc p~ 1 m ( - n+5 ^ 6 for the virial 
radius; u 2 ~ Gm/r v i T ~ /im' 1- "" 6 for the virial velocity; 
and p oc m/r 3 ir oc p 3 m~( n+3 ^ 2 for the average gas density. 
This suggests the average pressure scaling 

2 4 -2(n+2)/3 

p sc U OC p 771 " . 

We adopt n = -1.5, M = 1O 15 M , m = 3 x 10 13 M Q , 
and a bulk velocity of the subcluster of order the cluster 
velocity dispersion, v EC ~ lOOOkms -1 . For comparison, the 
merging components of the Coma cluster have virial masses 
0.9 x 1O 15 M and 6 x 1O 13 M (Colless & Dunn 1996). We 
also assume that the cluster and subcluster correspond to 
density fluctuations of the same value of p. Then Eq. J7J 
yields Ro = 2.3i? sc at the cluster centre, r = 0, and Ro = 
3AR SC at the cluster core radius, r = r c . So, according to 
this criterion, gas within 2-3 subcluster core radii will not 



be stripped as the subcluster falls, along a radial orbit, into 
a cluster which is about 30 times larger in mass. 

Further, we take the gas core radius to be proportional 
to the virial radius. Indeed, Sanderson and Ponman (2003) 
suggest that the gas core radius is about 0.1r v i r for clus- 
ters with temperature exceeding 1 keV, or the mass of a few 
times 1O 13 M0. Then the subcluster gas core radius is about 
(m/M)' n+5 ''' 6 ~ 0.13 times the cluster gas core radius for 
m/M = 0.03. For a rich cluster with the virial radius 3Mpc 
and the core radius ten times smaller, or r c = 300 kpc, we 
obtain the subcluster gas core radius as 7? sc ~ 40 kpc. This 
implies for such subclusters the stripping radius of at least 



Ro ~ 100 kpc . 



We adopt these values for qualitative estimates, keeping in 
mind that scatter about the fiducial values is likely to be 
large. 

Heinz et al. (2003) simulated a subcluster with a shal- 
lower gas density profile (f3 = 0.5), a somewhat large core 
radius R sc = 250 kpc, central gas number density p sc0 — 
3.6 x 10 -3 cm -3 and a temperature T = 3.2 keV, moving 
through a uniform gas of a density p c0 = 4.6 x 10~ 4 cm~ 3 
and temperature T — 7.7 keV. These authors find that the 
subcluster gas within Ro — 2i? sc survives ram pressure strip- 
ping, which compares favourably with our estimates based 
on Eq. JBJ. 

Flow past a solid sphere develops into a turbulent wake 
for sufficiently large Reynolds numbers. Experiments and 
numerical simulations (Tomboulides & Orszag 2000 and ref- 
erences therein) show that the transition to turbulence oc- 
curs at Re « 400, via the Kelvin-Helmholtz instability of a 
shear layer that results from the separation of the bound- 
ary layer on the sphere's surface. It is not clear what is the 
critical Reynolds number for a gaseous sphere. It can be 
speculated that the entrainment of the dense subcluster gas 
into the flow can be a cause of the flow randomness addi- 
tional to that past a solid sphere. The Kelvin-Helmholtz 
instability does indeed develop on the boundary between 
the subcluster and the ambient gas, e.g., in the simulations 
of Takizawa (2005) among many other authors, leading to 
prominent eddy-like structures in the subcluster wake that 
can be described as a turbulent flow. 

Nulsen (1982) describes how the introduction of eddies 
of a scale I can make the boundary layer smooth on this 
scale, suppressing the Kelvin-Helmholtz instability at wave- 
lengths smaller than I. Longer- wavelength modes are still 
unstable, and the largest unstable scales are comparable to 
the stripping radius. The Kelvin-Helmholtz instability is ef- 
ficient in eventually removing gas from the subcluster. Ac- 
cording to Heinz et al. (2003), all the gas is removed after a 
time of order a few times 10Ro/v c , that is a few billion years. 
This implies that a subcluster can generate a turbulent wake 
during one or two passages through the cluster. 

Altogether, the flow produced by the Kelvin-Helmholtz 
instability can produce turbulence in the subcluster's wake, 
subject to the same reservations as discussed after Eq. 
It is then plausible that the wake far downstream of the 
subcluster is well described by Prandtl's self-similar solution 
for turbulent wakes. 
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2.3.2 The area covering and volume filling factors of 
turbulent wakes 

Prandtl's solution for the turbulent scale and velocity vari- 
ation with distance x along the wake has the form (Landau 
& Lifshitz 1975) 



l x — Li(x/Li) 



1/3 



Vox ^ Vi(x/Li)' 



-2/3 



(8) 



where the turbulent velocity near the head of the wake can 
be identified with some fraction of the subcluster speed, 
Vi — v c — 1000 kms -1 ; and the initial value of the tur- 
bulent scale Li is close to the stripping radius Ro estimated 
above, Li ~ Ro. The Reynolds number varies along the wake 
as 

LiVi 



Re(x) = Rei(x/Ro) 



-1/3 



Rei 



and we assume that the wake remains turbulent as long 
as the Reynolds number exceeds a certain critical value, 
Re cr . Then the length of a turbulent wake, X, follows from 
Re(X) = Re cr as 



X_ ^ / Rei \ 
Ro ~ \R^~ T ) 



(9) 



where Re cr = 400 (Tomboulides & Orszag 2000) can be 
adopted for illustrative purposes. Thus, we assume that the 
critical value of Re required to maintain turbulence within 
the wake is the same as that to produce it immediately be- 
hind a solid sphere. This will result in a quite conservative 
estimate of the wake length since turbulence can plausibly 
sustain at even smaller local values of Re. 

The area of a single wake of a length X seen from the 
side is given by 

Ro{x/R,f ri dx=±Rl |^g-J -1 . 

If the wake axis is inclined by a random angle a to the line 
of sight, where a is uniformly distributed with < a < n, 
the average area of a single wake in the sky plane is given 
by 



sin 7 73 a da « 0.74 



The area covering factor of N wakes within a region larger 
in diameter 2r than the wake length, r > X/2, follows as 



fs 



NS 



(10) 



Similarly, the volume filling factor of N wakes is given by 
9 „ /i?o\ 3 17 X \ 5/3 



(11) 



Both estimates assume that the wakes do not intersect in 
three dimensions, which makes them slight overestimates. 

The r.m.s. turbulent velocity averaged over the wake 
length X is given by 



v 0x dx 



V3Vi 



o 

Ro 
X 



1/2 



Ro\ 1/3 



(§) 



1/2 



Re cr 
Re ; 



1/2 



(12) 



Similarly, the typical turbulent scale can be identified with 
the average wake width, 

3 



l ~ -R (X/Ro) 



1/3 



(13) 



The area covering factor and the volume filling factor of 
the wakes sensitively depend on the subcluster mass m and 
are larger for larger m (which is consistent with the idea 
that turbulence produced in major merger events fills the 
volume) . We also note the strong dependence of the covering 
and filling factors on the Reynolds number: fs oc Re 4 and 
fv oc Re 5 (for X/Ro > 1). 

It is clear from Eqs 1101 and li lt that one can have 
fv <C 1 but /s ~ 1 (note that Ro <C r). When this is 
the case, the wakes are widely separated in space, but their 
images overlap in the two-dimensional projection onto the 
sky plane. 



2.3.3 Subcluster wakes 

For subclusters of a mass m = 3 x 1O 13 M0, we adopt Ro = 
100 kpc, Vi = c s and A = 1 kpc to obtain, from Eq. © , an 
estimate 

3 / , , . \ -3 



X_ 
Ro~ 



27 



Ro 



100 kpc 



4XS 
1 kpc 



/RgcrV 

V 400 ) 



(14) 



With the scaling of Eq. {JSJ, the merger rate of subclusters of 
this mass is about 5 times larger than that of major mergers; 
thus, we assume that JV = 5 subclusters of this mass can 
(almost) simultaneously fall into a larger cluster. The area 
covering and volume filling factors of N — 5 wakes within 
the radius r = r v i r w 3 Mpc are estimated as 



/s~0.2£ 



N 



Ro 



100 kpc 



V 400 ) 



and 



0.02 



N f Ro 



5 V 100 kpc 



/Recr 
V 400 



4S\ 
1 kpc 



/ 4SX V 
^lkpcj 



The covering and filling factors strongly depend on the 
poorly known viscosity, parameterized with S. For S < 0.16, 
we obtain /s > 1, but the volume filling factor remains 
smaller than unity for <5>0.1. Furthermore, both fs and 
fv depend on high powers of another poorly known pa- 
rameter, the stripping radius Rq. Hence, properties of the 
subcluster wakes can be rather different in apparently simi- 
lar clusters. In addition, results of numerical simulations of 
turbulent wakes should be treated with caution as other- 
wise reasonable approximations, numerical resolution, and 
numerical viscosities can strongly affect the results. 

Upper limits on the covering and filling factors follow if 
we assume that the wake length is equal to or exceeds the 
region size, X — 2r v i r = 6 Mpc, which is obtained in Eq. JUJ 
if Re ; /Re cr > 4 (and A = 1 kpc): 



fs< 0.5(45)- 



fv < 0.08(45)" 



which yields f s < 2 and f v < 0.5 for 8 > 0.2. 

Thus, wakes from subclusters of a mass 3 x 1O 13 M0 can 
occupy just a small fraction of the total volume within the 
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viral radius of a cluster, but their area covering factor can 
be substantial. Given the sensitive dependence on the poorly 
known value of the Reynolds number, it appears reasonable 
to assume that fs = 0(1), that is any line of sight within 
the virial radius will have good chance to intersect at least 
one turbulent wake. 

The r.m.s. turbulent velocity averaged over the wake 
length follows from Eq. 112H as no — 260 km s" 1 if averaged 
along the whole length X ~ 2.7 Mpc, and vo — 190 km s -1 
within the cluster virial radius (with X = 2r v i r — 6 Mpc). 
The average turbulent scale follows from Eq. 113H as Iq ~ 
200 kpc. 

Takizawa (2005) has recently studied turbulence gen- 
erated by a subcluster of total mass of 1O 14 M , gas core 
radius of 100 kpc and a much higher central density ~ 
3 x 10 -2 cm -3 , moving through a uniform medium about 100 
times less dense. Turbulent velocities obtained in those sim- 
ulations, 300-500 km s _1 (see also Norman & Bryan 1999; 
Ricker & Sarazin 2001; Schueker et al. 2004), are in a rea- 
sonable agreement with our estimates. 

2.3.4 Galactic wakes 

The stripping radius of galaxies could be estimated similarly 
to that of subclusters, but the arguments are complicated 
by the replenishment of interstellar gas by stellar winds, 
magnetic fields that affect the Kelvin-Helmholtz instability, 
etc. Both numerical models (Portnoy et al. 1993; Balsara et 
al. 1994; Acreman et al. 2003) and observations (Sun et al. 
2005) indicate that gas within 

i?o=3-5kpc (15) 

of the centre of a massive elliptical galaxy can remain un- 
stripped. We assume that the Reynolds number based on 
this scale, Re ~ (10-15)5 _1 , is large enough to produce a 
turbulent wake, e.g., because S is small enough. 

Consider a rich galaxy cluster, where N ~ 100 galaxies 
are found within the gas core radius r c = 180 kpc (Sarazin 
1988). From Eq. 0, the 

area covering factor of galactic 
wakes in this region is unity if 

30-15, X~ 100-70 kpc, (16) 

Ho 

where the range corresponds to that in Eq. I|150 . Wakes of 
this length would require Rei/Re cr ~ 3, which is obtained, 
e.g., for S ~ 0.01 if Re cr = 400. 

The r.m.s. turbulent velocity and scale averaged along 
the wake follow from Eqs il 1 "2 fi and 1 Kit as 

v ~ 300kms _1 , Z ^8kpc (17) 

for Ro — 4 kpc and X/ Rq — 20. The volume filling factor of 
such wakes is fv — 0.07. 

The size of galactic wakes required to cover the pro- 
jected cluster area, given by Eq. 1161 . does not seem to be 
unrealistic. For example, Sakelliou et al. (2005) have ob- 
served a wake behind a massive elliptic galaxy (mass of or- 
der 2 x 1O 12 M ) moving through the intracluster gas at a 
speed about v c ~ 1000 kms -1 . The length of the detectable 
wake is about X ~ 130 kpc (assuming that it lies in the 
sky plane) , and its mean radius is 40 kpc (obtained from the 
quoted volume of about 2 x 10 6 kpc 3 ). These authors argue 
that the wake is produced by the ram pressure stripping of 



the interstellar gas. The projected area of the wake is about 
10 4 kpc 2 , as compared to 10 3 kpc 2 for the wake parameters 
derived above. This wake has been detected only because it 
is exceptionally strong, and it is not implausible that weaker 
but more numerous galactic wakes can cover the area of the 
central parts of galaxy clusters. 

We conclude that subcluster wakes are likely to be tur- 
bulent, but galactic wakes can be laminar if the viscosity of 
the intracluster gas is as large as Spitzer's value. Given the 
uncertainty of the physical nature (and hence, estimates) 
of the viscosity of the magnetized intracluster plasma, we 
suggest that turbulent galactic wakes remain a viable possi- 
bility. Both types of wake have low volume filling factor but 
can have an area covering factor of order unity. 



3 MAGNETIC FIELD IN THE 
INTRACLUSTER GAS 

In this section we discuss the amplification of an initially 
weak seed magnetic field by the fluctuation dynamo oper- 
ating in the intracluster gas. The seed field itself can be 
produced by a wide range of mechanisms (Appendix [X] see 
also Ruzmaikin et al. 1989; Widrow 2002; Brandenburg & 
Subramanian 2005). We first discuss the fluctuation dynamo 
in general terms. These general ideas are then applied to the 
various contexts of intracluster turbulence discussed above. 
First, we consider the merger epoch when the turbulence 
can be assumed to be in a statistically steady state, then 
the later epochs after the driving by the merger had ceased 
and the turbulence decays, and finally to magnetic field gen- 
eration in turbulent wakes. 

3.1 The fluctuation dynamo 

The evolution of a magnetic field embedded into a flow 
of conducting fluid is controlled by the magnetic Reynolds 
number defined, similarly to Eq. as 

_ vplp 

Am — ) 

n 

where 77 is the magnetic diffusivity (inversely proportional 
to the electric conductivity). 

The exponentially fast amplification of an initially weak 
magnetic field by a random flow (called the fluctuation dy- 
namo) is a result of a random stretching of magnetic field 
by the local velocity shear (see reviews in Zeldovich et al. 
1990 and Brandenburg & Subramanian 2005). For R m ^> 1, 
magnetic field is nearly frozen into the flow. Then, due to 
the random stretching, magnetic field lines grow longer, that 
is B/p increases, where p is the gas density. For flows with 
p approximately constant, the magnetic field will be am- 
plified. Such amplification comes at the cost of a decrease 
in the scale of field structures in the directions perpendicu- 
lar to the stretching (i.e., on average in all directions if the 
flow is statistically isotropic). This enhances Ohmic dissi- 
pation and the latter ensures that the correlation function 
of magnetic field can grow exponentially as an eigenfunc- 
tion if the Lorentz force is negligible (the kinematic dy- 
namo). The growth occurs under a fairly weak condition 
Rm > Rm,a — 30-100 (where the variation within the range 
depends on the form of the velocity correlation function). If 
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Vi is the velocity at a scale I, the e-folding time for the mag- 
netic field is roughly equal to the eddy turnover time I /vi . In 
the Kolmogorov turbulence, where vi oc l 1 ^ 3 , the e-folding 
time is shorter at smaller scales, l/vi oc l 2 ^ 3 , and so smaller 
eddies amplify the field faster. 

Since 77 <C v in the rarefied intracluster plasma (e.g., 
Brandenburg & Subramanian 2005), we have 7? m Re. 
Therefore, _R m 3> R m ,cr if Re > 100, so that random motions 
in galaxy clusters will be a dynamo for any Reynolds number 
which is large enough to make them turbulent. 

Numerical simulations of magnetic field evolution in 
turbulent flows confirm that the fluctuation dynamo ac- 
tion readily occurs in forced and convective turbulent flows 
(Meneguzzi, Frisch & Pouquet 1981; Cattaneo 1999; Hau- 
gen et al. 2003, 2004; Maron, Cowley & McWilliams 2004; 
Schekochihin et al. 2004), especially when _R m > Re. Such 
simulations are also able to follow the fluctuation dynamo 
into the non-linear regime where the Lorentz forces becomes 
strong enough to affect the flow as to saturate the growth 
of magnetic field. 

In the kinematic regime, the field is predicted to be 
intermittent, i.e., concentrated into structures whose size, 
in at least one dimension, is as small as the resistive scale 

l v = loRm /2 (18) 

in a single-scale flow (e.g., Ruzmaikin et al. 1989; Zeldovich 
et al. 1990). We emphasize that magnetic field at the small 
Ohmic diffusion scale is produced by the shear of the flow 
at a larger scale lo- 
in a turbulent flow, where a broad spectrum of motions 
is present, flow at each scale I would produce magnetic struc- 
tures at all scales down to the corresponding Ohmic scale. In 
the kinematic regime this would correspond to a set of eigen- 
functions, each with a distinct growth rate vi/l. The fastest 
growing eigenfunction is due to stretching by the smallest 
eddies with Rm(l) > Rm,a, where Rm{l) = R m {l/lo) 3/A - 
These are the viscous scale eddies, with I = l v = ZoRe -3 ^ 4 , 
provided i? m /Re > -R m , cr . However in the nonlinear regime, 
when the fastest growing mode saturates, larger scale modes 
could still grow. Since most of the kinetic energy is contained 
at the scale lo, the dominant magnetic scale could still be 
determined by dynamo action due to eddies of scale lo and, 
especially, by the subtle details of the dynamo saturation. 
We now discuss how the dynamo action could saturate. 

Nonlinear effects can modify the resulting magnetic 
structures, although it is as yet not clear in what way (cf. 
Haugen et al. 2003, 2004; Schekochihin et al. 2004). A sim- 
ple model of Subramanian (1999) suggests that the smallest 
scale of the magnetic structures will be renormalized in the 
saturated state to become 

l B =i loRnS , (19) 

instead of the resistive scale l v - This essentially happens via 
a renormalization of the effective magnetic diffusivity in the 
models of Subramanian (1999, 2003). In other words, it is 
suggested that the dynamo action can be saturated via a 
reduction of the effective magnetic Reynolds number down 
to its critical value for the dynamo action. . 

1 Saturation could also happen if the stretching properties of 
the flow are suppressed by the Lorentz force (Kim 1999). This 



Such results, however plausible they are, require further 
substantiation, e.g., by numerical simulations. Dynamo sim- 
ulations of Haugen et al. (2003, 2004) with v/rj = 1, where 
i? m ,cr ~ 35, show that magnetic energy per unit logarith- 
mic interval of k, kMk, has a broad maximum at k ~ 9, 
a scale about 6 times smaller than the forcing scale, but a 
factor of about four larger than the resistive scale given by 
k v = fco-Rm 2 ^ (20-30)fc for R m = 420-960, in agreement 
with the above idea and Eq. (1191 . (We note, however, that 
it is not quite clear how significant is the difference or agree- 
ment here since all the estimates have factors of order unity 
omitted, which can be important at the modest values of 
Rm available.) Further, these simulations also show that the 
value of fcs does not scale with i? m when i? m is increased 
from about 420 to 960, confirming Eq. 1191 . However, the 
magnetic spectrum is rather broad and it is difficult to iden- 
tify accurately the dominant magnetic scale in those simu- 
lations. Nevertheless, it is clear that the nonlinear magnetic 
field distribution is less intermittent (i.e., its scale is larger) 
than at the kinematic stage. Below we present evidence for 
this nonlinear behaviour in our own simulations of the fluc- 
tuation dynamo. The values of magnetic Reynolds number 
accessible now in such simulations are too modest to make 
any confident conclusions, but we believe that our approach 
to the saturation of the fluctuation dynamo is consistent 
with the evidence available. 

For the Kolmogorov turbulence, in the kinematic 
regime, the corresponding resistive scale for the marginal 
mode is predicted to be l v — loRm 3 ^ 4 (Subramanian 1997; 
Brandenburg & Subramanian 2005), when it is larger than 
the viscous cut-off scale (that is when l v > l„ = loKe~ 3 ^ 4 or 
Pm < 1). This scaling is different from that in Eqs 1181 and 
1191 . which apply to a single-scale flow, because the shear- 
ing rate now is (vo/lo)(l/lo)~ 2 ^ 3 at any scale I in the inertial 
range. For the marginal magnetic mode (neither growing nor 
decaying), the shearing is balanced by dissipation at I = l v 
which occurs at a rate rj/l^. In the intracluster gas, we gen- 
erally have < l v . However in the saturated state one may 
have Ib > lv for Re > R m ,cr- This would then suggest a 
different scaling, Ib — loRu?,lr instead of Eq 1191 . Never- 
theless, in both the simulations presented below and in real 
clusters the flow is not strongly turbulent, i.e., it has no 
extended Kolmogorov inertial range because the Reynolds 
number is not very large; so Eq. 1191 can remain a bet- 
ter approximation. There is some evidence for this from the 
simulations, in that Eq. 1191 agrees better with the wave 
number at which the magnetic spectrum peaks. We shall 
therefore use Eq. 1191 in our estimates. 

We note that properties of MHD turbulence can depend 
on the ratio P m = u/r), known as the magnetic Prandtl num- 
ber. The intracluster gas has P m 2> 1 if Spitzer's viscosity 
and resistivity are adopted; realistically large values of P m 
are not accessible to computer simulations, but we shall dis- 
cuss simulations with a modestly large value of P m in what 
follows. 



can perhaps be described as a reduction of the effective magnetic 
Reynolds number to its critical value i? m , C r of the kinematic dy- 
namo, a feature which however has not yet been studied in the 
model of Kim (1999) 
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3.2 Application to cluster turbulence 

Here we present semi-quantitative estimates to character- 
ize the fluctuation dynamo in the intracluster gas, before 
discussing, in Sect. |1J direct numerical simulations of the 
fluctuation dynamo. 



3.2.1 The epoch of cluster formation 

During the epoch of major mergers, we expect that the intr- 
acluster medium is involved in a steady-state, driven turbu- 
lence, for which we assume the Kolmogorov spectrum. Due 
to the action of eddies at scale I, the r.m.s. magnetic field 
grows exponentially at a rate 



7(0 



Hi = _L (JJ) 

l ~ toi UoJ 



-2/3 



tOi — , 



(20) 



where subscript 'i' refers to the initial, steady state of the in- 
tracluster turbulence. As summarized in the last paragraph 
of Sect. 12. ll the turbulent speed and scale can be adopted as 
voi — 300 km s -1 and Zoi = 150 kpc, respectively. Assuming 
that the driven turbulence lasts for U — 3 Gyr, we obtain an 
amplification exponent of the magnetic field T — j(l)t{ ~ 6 
at I — loi (and larger at smaller scales). So the seed field 
can be amplified by a factor 400 by motions at I = Zqi dur- 
ing this time (the amplification factor is larger at smaller 
scales). For a seed magnetic field of 10 -8 G, this amplifi- 
cation is sufficient to explain the observed magnetic fields; 
this implies that the observed magnetic fields are plausibly 
in the saturated state and the Lorentz force can now affect 
significantly the velocity field in galaxy clusters. 



3.2.2 The epoch of decaying turbulence 

After the driving forces have been diminished, the turbu- 
lence decays. Both the instantaneous outer scale and the 
viscous scale then increase with time (the latter, because of 
the decrease in the Reynolds number). In this situation it is 
more useful to estimate the growth rate at a fixed scale Z, 
which belongs to the turbulent spectrum for a long time and 
can hence lead to a fluctuation dynamo, instead of consider- 
ing an evolving viscous scale (where the dynamo time scale 
is the shortest). The instantaneous growth rate of the r.m.s. 
magnetic field due to motions at a fixed scale I decreases 
with time as 



'■'n 



(Z 2 Zo0 1/3 V toi 



t > t 0i 



(21) 



Here we have adopted a Kolmogorov spectrum with vi = 
Wo(t)[Z/io(t)] an d for numerical estimates take a = 6/5 
and P = 2/5. For a growth rate evolving with t, magnetic 
field evolves as B tx exp J y(l, t') dt'. If turbulence is main- 
tained in a steady state during an initial period t < tt and 
then 7 decreases as in Eq. 12 U . the amplification exponent 
for the Kolmogorov spectrum follows as 



10 

where 



7 (M')d*' = r(*)(£) 



-2/3 



(22) 



r(t) = 



t 

hi ' 



t<t[+ t oi , 
, t > t{ + toi , 



where £ = 1 — |a — |/3. We have applied the power-law 
12 1> only at t > U + toi and assumed that the decay does 
not affect the growth rate before that time. Assuming that 
the turbulence starts decaying after a time U = 3 Gyr, we 
obtain, at t = 5 Gyr, that the energy-range speed reduces 
down to vo ~ 130 kms , whereas the energy-range scale 
increases to lo « 260 kpc. The amplification exponent of 
magnetic field due to motions at I = loi is obtained as V ~ 9 
at t = 5 Gyr (consisting of 7 at t — tf + toi = 3.5 Gyr and 
only less than 2 at later times), so that the seed field could 
be amplified by a factor about 6 x 10 3 by the end of the 
decay phase if it were too weak to bring the dynamo to the 
saturated state earlier. In order to obtain magnetic field of 
1 /iG at this scale for t — 5 Gyr, a seed field of 2 x 10~ 10 G 
would be sufficient. Again the amplification is larger due 
to smaller scale eddies which remain part of the turbulent 
cascade as the turbulence decays. 

3.2.3 Dynamo action in wakes 

Using the turbulent speed and scale averaged over the wake 
length, as derived in Sect. 12.3.21 we obtain magnetic field 
growth time scales 7 -1 ~ lo/vo — 0.8 Gyr for subcluster 
wakes and 3 x 10 7 yr for galactic wakes. At a given posi- 
tion, time available for the dynamo action is X/Vi, where 
V\ ~ 1000 kms -1 is the speed of a subcluster or a galaxy. 
Therefore, the dynamo amplification exponent is given by 

r~^~3 

Vi lo ~ 

for both subclusters and galaxies, which implies additional 
amplification by a factor 20 at the outer scale. 



3.2.4 Magnetic field strength in the intracluster gas 

The maximum local magnetic field strength produced by the 
turbulent dynamo will be, presumably, close to equipartition 
with the turbulent energy: 



3a*G ( 



10- 



*f( 



200 kms 



(23) 



where p is the gas density. [We note that some models of 
nonlinear fluctuation dynamo predict stronger local mag- 
netic fields (Belyanin et al. 1993, 1994), but here we adopt 
a conservative limit 12311 .1 As discussed in Sect. 13.11 mag- 
netic field produced by the fluctuation dynamo is expected 
to be spatially intermittent (especially at early stages of 
dynamo action), i.e., represented by intense filaments and 
sheets whose volume filling factor is less than unity. Numer- 
ical simulations and analytical models recently reviewed by 
Brandenburg & Subramanian (2005) suggest that magnetic 
sheets and ribbons are prevalent, whose thickness is given 
by Eq. 119H and whose other two dimensions are of the order 
of the turbulent scale lo (see Fig.|SJl. Then the volume filling 
factor of magnetic structures within a single turbulent cell 
in the statistically steady state can be estimated as 
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fa — -^p — Rxa,lx w 0.17 , 
'o 

where the numerical value refers to P m ,cr = 35. Therefore, 
the r.m.s. magnetic field within a turbulent cell is of order 

{BY 2 - fTB„ 



100 



1.2 



MG (lO- 3 cm- 3 ) GoOkms" 1 ) 



TO" 3 cm" 3 

H n 1 



( Rm,cr \ 



(24) 



which implies that magnetic energy density in a saturated 

1/2 

dynamo state is a factor Rm,cx — 6 times smaller than the 
turbulent energy density. (Here and below angular brack- 
ets denote averaging.) We emphasize that weaker volume- 
filling magnetic fields are also present; their contribution to 
the magnetic energy density can be somewhat smaller than 
that of the intermittent part. In agreement with this esti- 
mate, magnetic energy density is about 0.2 (0.33) of the 
kinetic energy density in Model 1 (Model 2) of the numeri- 
cal simulations discussed in Sect.0 As discussed by Haugen 
et al. (2004, their Fig. 14), this ratio weakly varies with 
magnetic Reynolds number; in their simulations, the varia- 
tion is from about 0.25 for R m = 420 to 0.4 for P m = 960, 
both with P m = 1. For P m = 30, their simulations yield 
the magnetic-to-kinetic energy ratio of about unity. Thus, 
the flow at larger magnetic Prandtl number appears to be 
producing more magnetic energy; perhaps magnetic struc- 
tures fill turbulent cells more densely, and/or a smoothly 
distributed magnetic field is stronger. This feature is yet to 
be understood. What is, however, important for our imme- 
diate purpose here, is that the fluctuation dynamo produces 
magnetic fields whose energy density is comparable to the 
kinetic energy density of the turbulence. 

If the volume filling factor of the turbulent flow fy is 
less than unity, as in the case of turbulent wakes of sub- 
clusters and galaxies, the r.m.s. magnetic field in the cluster 
volume is obtained from Eq. (12 H by further multiplication 

1/2 

by a factor f v ; this is a measure of the total magnetic 
energy of the cluster. However, this quantity has little phys- 
ical significance because it would not result from any local 
magnetic measurement. In this sense, the local value 124H is 
more meaningful; it is presented in Table together with 
other quantities that characterize turbulence and magnetic 
fields at various stages of the cluster evolution. 



4 SIMULATIONS OF THE FLUCTUATION 
DYNAMO 

We have simulated the generation and subsequent decay 
of dynamo-active turbulence using the numerical model 
of the fluctuation dynamo by Haugen et al. (2003, 2004), 
where isothermal, viscous, electrically conducting, compress- 
ible gas is driven by a random force imposed as a source in 
the Navier-Stokes equation. The Navier-Stokes, continuity 
and induction equations are then solved in a Cartesian box of 
a size Dona cubic grid with 256 3 mesh points. The driving 
force / is sinusoidal in the spatial coordinates, transversal 
(/ _L k with k the wave vector of the force), and localized in 
the wave- number space about a certain wave number k — kf , 
so it drives almost incompressible, vortical motions in a cer- 
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Figure 2. The evolution of the r.m.s. fluid velocity (solid) and 
magnetic field (dashed) in driven and then decaying turbulence 
with k t = 5 (Model 1) (v rms = (v 2 ) 1 / 2 and likewise for B). 
Velocity is measured in the units of diffusive speed at the driving 
wavelength, ukf , and magnetic field is expressed in similar velocity 
units, (iirp) 1 / 2 c s vk{ . Hence, v Iuls numerically coincides with Re 
in the statistically steady state. Time is measured in the units 
of the initial turnover time of the energy-containing eddies, toi. 
Dotted line shows the asymptotics (t — tf) — °' 65 , where tf = is 
the time when the forcing is turned off. 



tain wavelength range around 2n/kf (see Haugen et al. 2004 
for details). The direction of the wave vector of the force 
and its phase change randomly every time step in the sim- 
ulations, so the force is effectively 5-correlated in time. 

We represent numerical results using the following 
units. (Tilde is used to denote dimensionless quantities.) For 
a unit length d, the computational domain size is equal to 
D — 2-ivd. The wave number is measured in the units of d~ . 
In simulations with dimensionless forcing wave number kf, 
it is appropriate to adopt kf = 2-k/Iq for its dimensional 
value, where lo = 150 kpc is the turbulent scale in a merg- 
ing cluster, as obtained in Sect. 12. ll Then the unit length is 
d — kflo I (27r) and the dimensional size of the computational 
domain is D — lok{. The unit density po can be adopted to 
correspond to the number density of no = 10 -3 cm -3 . The 
unit speed is the speed of sound, c s = lOOOkms -1 , so that 
the unit magnetic field is (Anpo) 1 ^ 2 ^ = 15 fiG. 

Here we report results obtained with two values of 
the central driving wave number kf. Some results were ob- 
tained with driving covering the range of dimensionless wave 
numbers k = 4.5-5.5, centred at kf = 5 (Model 1). Re- 
sults at higher resolution (which was especially needed when 
P m > 1), were obtained with the driving wave- number range 
of k = 1-2 centred at kf = 1.5 (Model 2). In these latter 
runs, the computational box contains just a few turbulent 
cells. 

The intensity of the driving was adjusted to obtain the 
r.m.s. Mach number of the turbulence of about 0.1 which 
produces relative density fluctuations of order 0.01 (imply- 
ing that only a small fraction of the total velocity is com- 
pressible). The kinematic viscosity and magnetic diffusivity 
in most runs are adopted to be equal to v = r\ — 2 x 10~ 4 c s d 
(i.e., a magnetic Prandtl number of unity). This corresponds 
to Re = Rm w 110 in Model 1 and Re = P m w 420 in 
Model 2, which is close to what is expected for Re in the 
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Table 1. Summary of turbulence and magnetic field parameters at various stages of cluster evolution: 
duration of the stage (the last two stages represent steady states), the r.m.s. velocity vo and scale ?o 
of turbulence and eddy turnover time to (for the decaying turbulence, values for the middle of the 
decay stage are given, 2 Gyr after its start), the equipartition magnetic field _B cq given by Eq. 1231 
(i.e., maximum field strength within a turbulent cell), thickness of magnetic filaments and sheets lg, 
defined in Eq. 1191 . for the statistically steady state of the dynamo, the r.m.s. magnetic field within a 
turbulent cell B rms = (B 2 ) 1 / 2 given by Eq. 1241 (the latter two obtained for -R m , C r = 35), and finally the 
standard deviation of the Faraday rotation measure ctrm [calculated using Eq. 1281 for the volume filling 
turbulence within 500 kpc of the centre and path length of 750 kpc in the first two lines, and assuming 
one transverse wake along the line of sight in the last two lines, using Eq. ESU- Subcluster mass of 
3 X 10 13 Mq has been assumed. 
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\pG] 


[kpc] 


[MG] 


[rad m~ 


Major mergers 


4 


300 


150 
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260 


2.0 


2 


44 
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Galactic wakes 




300 
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0.03 
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1.4 


1.6 


5 



intracluster gas. (We note, however, that the values of mag- 
netic Reynolds number explored here still are smaller than 
those expected in reality.) We have also considered the case 
where magnetic diffusivity is 30 times smaller than kine- 
matic viscosity in runs with hi = 1.5, v — 1.5 x 10 _3 c a d and 
rj = z//30, i.e., P m = 30, Re « 44 and R m w 1300.. Results 
presented in what follows refer to the case P m = 1 unless 
stated otherwise. 

In order to simulate dynamo action in forced and 
then decaying turbulence, the flow had been driven until 
it reached a statistically steady state, with a weak magnetic 
field introduced at the start of the simulation. Then the sys- 
tem was evolved for some period (about 45 time units in 
Model 1), after which the driving force was switched off; 
tf = is the time when the driving halts. The initial, weak 
magnetic field is random, with energy density of about 0.6% 
of the kinetic energy density in Model 1. 

Figure |5] shows the time evolution of the suitably 
normalized r.m.s. velocity and magnetic field obtained 
in Model 1, after the driving was switched off [in fact, 
v T m S /(vk{), shown with solid line, is the Reynolds number 
based on the forcing scale, with vo ~ frms]- The initial ex- 
ponential growth of the r.m.s. magnetic field that obtains 
(not shown in Figure |2j, is followed by its saturation at a 
level where its energy density is about 0.2 (0.3) of the tur- 
bulent energy density in Model 1 (Model 2). More precisely, 
the r.m.s. values of the turbulent velocity and magnetic field 
(measured in velocity units) in the steady state of Fig.|5|are 
about 0.114c s and 0.050c s , respectively, whereas the similar 
quantities for k{ — 1.5 are 0.116c s and 0.065c s . The critical 
value of the magnetic Reynolds number remains about 35 
in both models. 

The subsequent decay of both the velocity and magnetic 
field strength can be approximated by (t — if)" 0,65 for t 3> 
toi, as shown with dotted line. This decay law is consistent 
with Eq. |0J which predicts the power law exponent of —3/5. 
However, the alternative value a — 10/7 would result in the 
exponent of —5/7 ~ —0.71, which is also consistent with our 
numerical results. 

With lo — 150 kpc and n — 10~ 3 cm -3 , the r.m.s. turbu- 
lent velocity and magnetic field strength in the steady state 
in Fig. Hare {v 2 ) 1/2 « HOkms" 1 and {B 2 } 1/2 « 0.7/iG, 
respectively. These results favourably agree with estimates 



presented in Tabled in the sense that in both cases the ra- 
tio of the turbulent and AlfVen speeds is about 1/2, which 
confirms our estimate of the r.m.s. magnetic field strength 
in Eq. 12411 . In other words, if our simulations had stronger 
driving to achieve (u 2 ) 1//2 « 300 kms" 1 , then the r.m.s. mag- 
netic field would be (B 2 ) 1/2 « 2 /^G, as in our analytical 
estimate. 

The magnetic and kinetic energy spectra are shown in 
Fig- El I n the statistically steady state (the upper curves), 
kinetic energy in Model 1 peaks at fcf = 5, the driving wave 
number. However, magnetic energy has broad maximum at 
a significantly smaller scale, apparently because of its in- 
termittent structure. This difference is better visible in the 
right panel that refers to Model 2 where we have higher res- 
olution. Similar simulations, but with a significantly higher 
resolution 1024 3 (Haugen et al. 2003, 2004), confirm that the 
magnetic energy per unit logarithmic interval in the fc-space, 
feMfe, has a maximum at ks ~ 6fcf in excellent agreement 
with Eq. 1191 with i? m , cr = 35, but no significant depen- 
dence of kg on Rm has been revealed. The value of Rm is 
those simulations is about 1000, so the deviation from the 

— 1/2 

scaling ks ~ k v oc R m advocated by Schekochihin et al. 
(2004) is by a factor of 3-5 over the range R m = 400-1000. 

The length scales of the velocity and magnetic fields 
can be characterized more precisely in terms of their integral 
scales 



L v = 2tt 



Jk^Ek dk 
J E k dk 



(25) 



and similarly for the magnetic scale Lb; here integration 
extends over the whole interval of k available. These scales 
are simply related to the longitudinal Zl and transverse Zn 
integral scales of the magnetic fields by Zl = = \Lb 
[Eq. (12.91) of Monin & Yaglom 1975], and similarly for 
L v , but only approximately because v is not solenoidal. The 
dimensionless value of the longitudinal integral scale of mag- 
netic field in the steady state is then obtained from Fig.^Jas 
l B = w 0.16 (see AppendixIQ. This agrees reason- 

ably well with the prediction from our heuristic estimates, 
Ib ~ (2Tv/k F )Rm% 2 ~ 0.2. The time variation of these scales 
is shown in Fig. 2] When the turbulence decays, the integral 
scales of both velocity and magnetic field exhibit power-law 
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Figure 3. The energy spectra of kinetic (solid) and magnetic (dashed) energies at various stages of evolution, t/toi = 0, 10, 50, with 
(a) fef = 5 — Model 1, left panel, and (b) fcf = 1.5 - Model 2, right panel. Spectra obtained at later times are at lower levels because of 
the decay of turbulent energy. 
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Figure 4. Evolution of the integral scales of the velocity (L„, 
solid line) and magnetic fields (Lb, dashed line), as defined in 
Eq. with fc f = 5 (Model 1). 



increase, in agreement with Eq. (|1J; the growth slows down 
when L v has grown to become comparable with the box size. 

Magnetic energy at small scales has, at early times, ex- 
cess over kinetic energy because magnetic field is very in- 
termittent, which is especially clearly visible in the right 
panel of Fig. At later stages, magnetic field distribution 
becomes more homogeneous and this feature disappears. Si- 
multaneously, the scale of magnetic field increases and be- 
comes comparable to that of the flow, which is not the case 
at early stages. 

Figure illustrates (using Model 2) the structure of 
magnetic field in a turbulent flow in a statistically steady 
state (left panel) and at a late stage of decay (right panel). 
The magnetic field produced by the fluctuation dynamo 
consists of an intermittent part, represented by randomly 
distributed, intense magnetic ribbons, sheets and filaments 
(which can even be folded), immersed in a sea of volume- 
filling random magnetic field. The intermittency gradually 
reduces as the turbulence decays together with magnetic 
field because structures of smaller scale decay faster, and the 



volume filling factor of magnetic field increases with time - 
this tendency can easily be seen in the right panel of Fig. [S] 



5 OBSERVATIONAL DIAGNOSTICS 

5.1 The Faraday rotation measure 

An important observational diagnostic of the intracluster 
magnetic field is the Faraday rotation of polarized radio 
emission of background sources (located beyond the clus- 
ter or in its centre) produced in the intracluster gas. The 
Faraday rotation is quantified by the Faraday rotation mea- 
sure 

RM = K J n c B ■ dl , (26) 

where n c is the number density of free thermal electrons, the 
integral is taken along the path length L from the source 
to the observer, and K = 0.81 radm~ 2 cm 3 fiG" 1 pc _1 . For 
a magnetic field with zero mean value (B) = 0, the mean 
value of RM vanishes, whereas its standard deviation can be 
represented in the form (Appendix iBl see also Burn 1966; 
Sokoloffet al. 1998), 



CRM 



RMovJV 



(27) 



where RMo is the Faraday rotation measure produced in a 
single turbulent cell of a size lo and N = L/lo is the number 
of the cells along the line of sight. 

Suppose that each turbulent cell contains randomly ori- 
ented magnetic sheets of thickness Ib where magnetic field 
strength is equal to B eq , with a covering factor of order 
unity, as described in Sect. 13.2.41 Then the Faraday rota- 
tion measure produced in a single turbulent cell follows as 



RMo — Kn e B e qlB, and, adopting Ib = loRi 
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(28) 
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Figure 5. Snapshots of magnetic field in a cross-section through the middle of the computational domain of Model 2 at t/toi = 0.30 
when the system is in a statistically steady state (left panel), and t/toi = 59.34, at a late stage of decay (right panel). Here fcf = 1.5, so 
each frame contains a few turbulent cells. The magnitude of the field component perpendicular to the plane of the figure is shown color 
coded (in shades of grey) with black corresponding to field pointing into the figure plane, and lighter shades, to field pointing out of the 
plane. The field in the plane of the figure is shown with vectors whose length is proportional to the field strength. 



If the magnetic sheets multiply cover the projected area a 
turbulent cell, by a factor q = O(l), then additional fac- 
tor q 1 / 2 has to be included in Eq. 1281 . The comparison of 
Eq. 12811 with both numerical simulations and observations 
of Faraday rotation in galaxy clusters suggest that q m 1. If 
B cq oc vq oc t~ 3 ^ 5 and lo oc t 2 ^ 5 , the observed RM will, on 
average, decrease with time as 

o-rm oc [(t-tt)/t 0i }- 2/5 . (29) 

We have calculated the Faraday rotation measure for 
256 2 lines of sight through our computational domain, and 
then computed ctrm as the standard deviation of the results. 
The evolution of ctrm in decaying turbulence is shown in 
Fig. HJ for Model 1 and it exhibits remarkable agreement 
with Eq. 1291 . At earlier stages of the simulations in Model 1, 
crm first grows rapidly while magnetic field is exponentially 
amplified, and then remains fairly constant, <trm ~ 0.3, for 
—20 < t < 1 (Fig. [SJ. In dimensional units, this corresponds 
to ctrm ~ 80radm -2 . Model 2 results in a value of ctrm ~ 
0.47 (Fig. |HJ| , corresponding to 130radm -2 . 

It is useful to compare results of the simula- 
tions with the analytical estimate of Eq. 1281 . ctrm — 
Ly 2 in dimensionless units used in 
Fig. HO This gives ctrm ^ 0.4 for fc f = 1.5 in Model 2 and 
ctrm — 0.75 for kf — 5 in Model 1. This estimate of ctrm for 
Model 2, which has a higher spatial resolution, is in good 
agreement with the numerical simulations, but that obtained 
for Model 1 is a factor of about 2 lower than expected. Nev- 
ertheless, our simulations confirm that Eq. 128H and Tabled 
provide reasonably good estimates of the expected amount 
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Figure 6. The width of the histogram of Faraday rotation mea- 
sures, ctrm; calculated along 256 2 lines of sight through the com- 
putational box, with fcf = 5 (Model 1, as in Fig.|5J, as a function 
of time. The solid line is the least square fit to the data points 
at t > 1, while the dotted line corresponds to t -2 / 5 . Here Ctrm 
is measured in the units Kn e (B 2 ) 1 ^ 2 lo 280radm — 2 , so that 
ctrm ~ 80radm — 2 in the steady state. 



of Faraday rotation by magnetic field generated by the fluc- 
tuation dynamo. 

Altogether the estimate l|28|l and the amount of Faraday 
rotation in our simulations agree very well with observations 
of Faraday rotation in the intracluster gas. 
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Figure 7. The normalized autocorrelation function c(R) = 
C( R )/°RM> where C(R) = (RM{X + R)RM(X)) and cr 2 M = 
(RM 2 (X)) obtained for Model 2 (k { = 1.5) with P m = 1, at var- 
ious times: t = (solid), t = 30 (dashed) and t = 70 (dotted). 
The former curve refers to the statistically steady state, whereas 
the latter two illustrate how RM distribution becomes less inter- 
mittent as turbulence decays. Here R is measured in the units of 
k f l /(2w) Kb 70kpc. 
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Figure 8. The histogram (probability density) of the Faraday 
rotation measure calculated along 256 2 lines of sight through the 
computational box, with k{ = 1.5, for P m = 1 (solid), P m = 30 
(dashed) and P m = 1/4 (dash-dotted), all in the statistically 
steady state before the decay starts, and with Re = 450, 44, 445 
respectively. Here RM is normalized by Kn c {B 2 ) 1/l2 lo, with 
Brms = (P 2 ) 1//2 . 



FigureQshows the autocorrelation function of the Fara- 
day rotation measure for Model 2 at the beginning of the 
evolution and at two later times. As turbulence decays and 
magnetic field becomes less structured, the correlation scale 
of RM fluctuations increases. In the steady state, the di- 
mensional value of the Taylor microscale (or differential 
length scale) of the RM fluctuations is, as expected, about 
-Rrm ~ 15kpc, i.e., about half the thickness of magnetic 
sheets, Ib, as quoted in Table Q 

The situation is somewhat different for the wakes of sub- 
clusters and individual galaxies. As discussed in Sects |2~!01 
and 12.3.41 their volume filling factor is small, whereas the 
area covering factor can be of order unity. In other words, a 
line of sight typically passes through just a single turbulent 
wake, where the turbulent scale is comparable to the wake 
width. The resulting Faraday rotation measure is given by 

ctrm = Kn c B aq lB 

~ 8 - 6 nS ( 10-3^-0 (5g) (l^fc) ' (30) 

where we retain the notation ctrm because the resulting 
amount of Faraday rotation will remain random, both be- 
cause of the random orientation of the wakes and due to the 
randomness of magnetic field within the wake. 

Estimates of typical Faraday rotation measures ob- 
tained from Eqs 1281 and 1301 are given in the last column 
of Table □ 

Schekochikhin et al. (2005a, b) suggest a different spatial 
structure of the cluster magnetic field, based on their inter- 
pretation of the fluctuation dynamo for P m 3> 1 (Schekochi- 
hin et al. 2004). These authors suggest that magnetic field 
produced by the dynamo is locally anisotropic and repre- 
sents Zo-long magnetic sheets and/or ribbons multiply folded 
at a microscopic resistive scale l v . Such a model would typ- 
ically produce a much smaller ctrm- Arguments similar to 
those that lead to Eq. 1281 . (even assuming that the folds are 



randomised by some unspecified mechanism) predict ctrm 
smaller than in our model by a factor F ~ (P m / Rm,^) 1 ^ 2 ■ In 
our simulations, this corresponds to F ~ 3-5 over the range 
Rm = 400-1000. So the simulations that we have analysed 
(with Rm ~ 400), cannot help to confidently discriminate 
between the two magnetic field geometries. 

However, observations indicate that magnetic coherence 
scale is at least a few kpc and more plausibly exceed 10 kpc. 
This would be difficult to produce in the model of Schekochi- 
hin et al. (2005a,b), unless the effective value of -R m in the 
clusters is reduced, say due to plasma effects, to be close to 
Rm.ci- Furthermore, Schekochihin et al. (2005a,b) envisage 
systematic reversals of the folded magnetic field along the 
line of sight, rather than random changes of its direction. 
Such a systematic behaviour would reduce ctrm even fur- 
ther due to systematic cancellations of magnetic field along 
the line of sight which would preclude the random walk of 
the polarization angle assumed in Eq. I28H . 

We have also used numerical simulations to examine the 
effects of varying the magnetic Prandtl number on magnetic 
field structure and Faraday rotation. The probability distri- 
bution of the Faraday rotation measure along 256 2 lines of 
sight through the computational box is shown in Fig. |^1 for 
three values of the magnetic Prandtl number. The shape of 
the probability distribution is close to a Gaussian curve for 
P m = 1/4 and P m = 1 (which is a parabolic shape in this 
representation), but the distribution obtained at P m = 30 
exhibits shorter tails at large |RM|. The reason for this is 
apparently the abundance of small-scale structures that pro- 
duce smaller Faraday rotation when P m > 1, i.e., when the 
magnetic dissipation scale is smaller. Nevertheless, the stan- 
dard deviation of the Faraday rotation measure has similar 
values for both P m = 1 and P m = 30, ctrm ~ 0.47 and 0.3 in 
the units of Fig. |H1 respectively. This implies that magnetic 
field does not become more strongly folded as P m increases. 
For the reader's convenience, we note again that Re w 44 
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Figure 9. As in Fig.|3 but for different magnetic Prandtl num- 
bers; P m = 1 (solid), P m = 30 (dashed) and P m = 1/4 (dotted), 
all in the statistically steady state with k{ = 1.5 before the decay 
starts. 



and R m » 1300 for P m = 30 and Re » 445 and R m rs 111 
for P m = 1/4. 

Figure[§]shows the autocorrelation function of the Fara- 
day rotation measure for various values of P m . The correla- 
tion scales for both P m = 30 and P m = 1 are comparable, 
but that which obtains when P m = 1/4 is a factor of two 
larger. 

We emphasize again that both the form of the RM 
probability distribution and the correlation function do not 
change much as P m increases from unity to 30. This suggests 
that our results can be robust and directly comparable with 
observations even though P m ^> 1 in the intracluster gas. 
Nevertheless, it would be important to clarify this issue fur- 
ther using simulations with higher resolution. Note that the 
Faraday rotation measure is proportional to the product of 
magnetic field and its scale; therefore, it remains dominated 
by large scales for magnetic spectra Mk oc k K with k < 1. 

5.2 Polarization of cluster radio halos 

Estimates of the scale of magnetic field in galaxy clusters 
obtained above are somewhat larger than what is usually 
adopted. We predict that the number of turbulent corre- 
lation cells along a path length of 750 kpc can be only 
about 3-5. As discussed above, magnetic field in each cell 
has an intermittent component of randomly oriented sheets 
and somewhat weaker fields in a volume filling component 
(see Fig. It can be expected that synchrotron emission 
produced in such a random magnetic field will be signifi- 
cantly polarized (assuming that magnetic field is well or- 
dered within individual magnetic sheets). As discussed by 
Sokoloff et al. (1998, their Sect. 5.1), the expected degree of 
polarization is a random quantity whose standard deviation 
is about 

a p ~ p N~ 1/2 , 

where po ~ 70% is the intrinsic degree of polarization and 
N is the number of turbulent cells within the beam cylinder, 
which yields a v ~ 30%, neglecting beam depolarization (see 
below). (We note that this estimate strictly applies if N 3> 



1.) The polarized emission would be confined to elongated 
structures (cross-sections of magnetic sheets) of Is = 20- 
40 kpc in width and Iq = 150-300 kpc in length. The intrinsic 
polarization plane should be perpendicular to the major axes 
of the elongated synchrotron structures since magnetic field 
is mostly parallel to the magnetic sheets. 

However, the fractional polarization observed from clus- 
ter radio halos is less than 2-10% at the wavelength A = 
21 cm (L. Feretti 2003, private communication; Govoni & 
Feretti 2004). No significant diffuse polarized emission in 
the Coma cluster has been detected by Thierbach, Klein 
& Wielebinski (2003) at wavelengths All. 2 cm and A6.2cm. 
The depolarization can be attributed to internal Faraday 
dispersion by the random magnetic field, where the degree 
of polarization will be further reduced to [Eq. (34) of Sokoloff 
et al. 1998] 

l-exp(-S) 4 2 

p = o P , S = 2A o-rm 

(which is strictly applicable when N 3> 1). Faraday disper- 
sion readily explains the lack of polarization at A21 cm where 
this equation yields p w 0.2% for <trm = 200radm -2 . The 
Faraday depolarization is weaker at shorter wavelengths, 
with p w 3% at All cm and p « 20% at A = 6 cm. How- 
ever, the linear resolution of the observations of Thierbach 
et al. (2003) was W = 110 kpc at All. 2 cm and W = 60 kpc 
at A6.2cm. Given that the thickness of the elongated polar- 
ized structures is of order Ib = 25 kpc, beam depolarization 
would further reduce the degree of polarization at least by a 
factor W/Ib to 0.5% at All. 2 cm and 8% at A6 cm. These es- 
timates indicate that the polarization of cluster synchrotron 
halos should be weak but detectable at sufficiently high reso- 
lution and short wavelengths. In reality, each correlation cell 
may contain a few magnetic sheets with independent direc- 
tions of magnetic field (cf. Fig. |KJ. Therefore, the effective 
number of magnetic sheets along the path length (and/or 
within the telescope beam) can be a factor 2-3 larger than 
adopted above and our values of the degree of polarization 
can be overestimated by a factor of two. Further polariza- 
tion observations of cluster radio halos at short wavelengths 
can reveal magnetic structures suggested here. 

Shear in the gas motions at a scale of a few hundred kilo- 
parsecs, produced during major merger events, can make the 
random magnetic field locally anisotropic. The anisotropy 
can also be produced by differential rotation and/or inho- 
mogeneous inflow in cluster cores. Anisotropic random mag- 
netic field can produce significant polarization of the syn- 
chrotron emission, with the polarization vector orthogonal 
to the direction of the maximum r.m.s. field strength (Laing 
1981; Sokoloff et al. 1998). This polarization can be observ- 
able if the shear regions are large enough as to avoid the 
cancellation of polarization along the line of sight or across 
the telescope beam. 

Govoni et al. (2005) report detection of polarized emis- 
sion from filamentary structures in the cluster A2255, of 
a size 200 kpc x 500 kpc (see also Murgia et al. 2004) , but 
the orientation of the polarization plane mostly disagrees 
with the above suggestions, unless the amount of foreground 
Faraday rotation is larger than that assumed by Govoni et 
al. 

Another situation where significant polarization of syn- 
chrotron emission in the cluster environment can be ex- 



Evolving turbulence and magnetic fields in galaxy clusters 17 



pected are the wakes of subclusters and galaxies. Since in- 
dividual lines of sight pass through one (or a few) wakes 
wherein the turbulent scale is comparable to the wake width, 
polarization due to the random magnetic field can be de- 
tectable. 



6 DISCUSSION 

There is growing direct and indirect evidence for the pres- 
ence of random - and plausibly turbulent - motions in the in- 
tergalactic gas of galaxy clusters. We have identified several 
stages in their evolution, from a statistically quasi-steady 
motion during the epoch of major mergers, to the stage 
of decaying turbulence that follows, and to a state where 
turbulence is confined to the wakes of relatively small sub- 
clusters and individual galaxies. Typical parameters of the 
velocity and magnetic fields at various stages of the cluster 
evolution are summarized in Table random velocities of 
«o = 150-300 kms -1 can be maintained at various stages of 
the evolution, and their scale is expected to be lo = 150- 
300 kpc, with the exception of galactic wakes where it can 
be of the order of 10 kpc. 

It is not quite clear whether or not the random motions 
in the intracluster gas can evolve into developed turbulence. 
This depends on the value of the Reynolds number, a mea- 
sure of the relative strength of nonlinear hydrodynamic ef- 
fects and, therefore, of the strength of the spectral energy 
cascade. If the flow remains laminar, the motions can decay 
faster after the end of the major mergers, and the wakes 
can have properties different from those discussed above. 
However, our numerical simulations suggest that the power- 
law decay establishes itself even for the Reynolds number as 
modest as Re ~ 100. 

The turbulent flow of magnetized gas can accelerate rel- 
ativistic particles required to produce cluster radio halos 
(Tribble 1993b; Brunetti et al. 2004; Cassano & Brunetti 
2005). Turbulent mixing in clusters has also been invoked 
in modeling the transport of heat (Cho et al. 2003; Kim & 
Narayan 2003; Voigt & Fabian 2004; Dennis & Chandran 
2005), and metals (Rebusco et al. 2005). The dissipation of 
the turbulent energy can help to balance the cooling of clus- 
ter cores (Fujita, Matsumoto & Wada 2004; Rebusco et al. 
2005). Clearly, turbulence in clusters seems to be useful to 
understand several diverse aspects of cluster physics. 

Random motions during and immediately after the ma- 
jor merger epoch are plausibly volume-filling. However, ran- 
dom flows produced by the wakes can have area covering fac- 
tor of order unity, but the volume filling factor of such wakes 
can be relatively small, fv — 0.1, leaving large quiescent re- 
gions between the wakes. Therefore, a typical line of sight 
passes through a turbulent region with an r.m.s. speed of a 
200-300 kms" 1 , producing some observational signatures of 
developed turbulence, and yet there is enough space to ac- 
commodate well-ordered morphological features apparently 
unaffected by any random motions. For example, in the core 
of the Perseus cluster, where the density is higher (Chura- 
zov et al. 2004), with n c > 10~ 2 cm" 3 for r < 100 kpc, the 
mean free path is smaller A < 0.5 kpc, and the wakes can 
have a potentially larger filling factors fv < 0.5, for 8 — 0.2. 
Thus the presence of long Ha filaments observed by Fabian 
et al. (2003, 2005) in the core of the Perseus cluster may 



not be inconsistent with various evidence for random mo- 
tions in this cluster core (Churazov et al. 2004; Rebusco et 
al. 2005). A possible signature of such spatially intermittent 
turbulence could be a specific shape of spectral lines, with 
a narrow core, produced in quiescent regions, accompanied 
by nonthermally broadened wings. It would be interesting 
to pursue this idea further in a more quantitative fashion. 

In the presence of random motions, any pre-existing 
magnetic field will be rapidly destroyed owing to a reduc- 
tion of its scale by the velocity shear. Even in a quiescent 
medium, any nonuniform magnetic would decay by driv- 
ing motions whose kinetic energy can be efficiently con- 
verted into heat because the intracluster gas is expected to 
be rather viscous. Therefore, random magnetic fields confi- 
dently revealed in many clusters through their Faraday ro- 
tation must be constantly maintained even if the electrical 
conductivity of the intracluster gas is large. 

However, the same random motions - either turbulent 
or not - will generate magnetic fields via the fluctuation dy- 
namo action at all the stages of the cluster evolution. The 
field is amplified by random shear, which reduces its scale 
along the directions perpendicular to the shear layers. This 
makes the spatial distribution of the magnetic field inter- 
mittent. Numerical simulations give an impression of strong 
field regions being largely confined into magnetic sheets and 
ribbons (and, with lower probability, filaments) wherein its 
strength is similar to that given by energy equipartition with 
the overall kinetic energy density. Following Subramanian 
(1999), we argue that the volume filling factor of the mag- 
netic structures within a turbulent cell (provided they are 
mostly sheets rather than filaments) is of order 0.1-0.2. Our 
numerical simulations confirm this picture, but add to it a 
weaker volume-filling magnetic background, so that the to- 
tal magnetic energy density is about 1/5-1/3 of the kinetic 
energy density of the random flow. The limited experience 
available with fluctuation dynamo models at large magnetic 
Prandtl number seems to indicate that magnetic fields can 
be closer to energy equipartition with turbulence in more 
realistic models. 

The (random) Faraday rotation measures produced by 
such magnetic fields are in the la range of 100-200 rad m -2 
in agreement with observations. We note, however, that, ac- 
cording to our estimates, the scale of the magnetic field is 
Ib = 20-40 kpc, i.e., a factor of a few larger than what is 
usually assumed. The scale of the field is smaller, of order a 
kpc, for galactic wakes; it could also be smaller if there were 
other sources of stirring like radio galaxies. The maximum 
field strength in the magnetic structures is about 2-4 ^iG, 
whereas its r.m.s. value within a turbulent cell is 1-2 fj,G. 
Such r.m.s. field strengths are in better accord with those 
inferred from synchrotron intensity assuming equipartition 
between magnetic fields and cosmic rays, or with inverse 
Compton limits. 

We predict that synchrotron emission from cluster radio 
halos similar to that in the Coma cluster can be significantly 
polarized at short wavelengths A = 3-5 cm. 
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APPENDIX A: SEED MAGNETIC FIELDS IN 
GALAXY CLUSTERS 

There is a number of sources of seed magnetic fields in galaxy 
clusters. It is well known that the intracluster medium 
(ICM) has high metallicity which must have been produced 
in stars in galaxies and subsequently ejected into the galactic 
interstellar medium (ISM) and then into the ICM. Since the 
ISM is likely to be magnetized with fields of order a few fiG, 
this would lead to a seed field in the ICM. One can roughly 
estimate the seed field resulting from stripping the galactic 
gas, by using magnetic flux conservation under spherically 
symmetric expansion; that is B ae ed — (piCM/piSM) 2 ^ 3 -B ga i. 
For B ga i ~ 3fiG, and picm/pism — 10~ 2 -10~ 3 , one gets 
S seB d ~ 0.1-0.03 pG. One may get even larger seed fields if 
there is a substantial number of active galaxies with magne- 
tized outflows: if about 10 3 galaxies have mass outflow with 
M ~ O.lMQyr -1 lasting for 1 Gyr, with a Poynting flux 
about 10% of the material flux, and the field gets mixed 
into the cluster gas over a Mpc sized region, B scc d — 0.3 [iG 
would result (Brandenburg 2000). This estimate, however, 
assumes that all the intracluster gas has been processed 
through the outflows, which may be an exaggeration. 

Another source of seed fields is likely to be the outflows 
from earlier generation of active galaxies (radio galaxies and 
quasars) (Rees 1994; Goldshmidt & Rephaeli 1994; Medina- 
Tanco & Enfilin 2001; Furlanetto & Loeb 2001; Colgate, Li & 
Pariev 2001). Such outflows can produce magnetized plasma 



bubbles in some fraction of the intergalactic volume (typi- 
cally of order 10% - Furlanetto & Loeb 2001) which, when 
incorporated into the ICM, would seed the general cluster 
gas with magnetic fields. If one assumes that the cluster 
gas is 10 3 times denser than the intergalactic medium and 
blindly uses the enhancement of the bubble field due to com- 
pressions during cluster formation, one can get fields as large 
as 0.1-1 fiG in the ICM (Furlanetto & Loeb 2001). However 
this is to ignore the issue of how the field in the magnetized 
bubble, especially if it is predominantly relativistic plasma 
from a radio galaxy, mixes with the unmagnetized and pre- 
dominantly thermal gas during cluster formation, and the 
resulting effects on both the field strength and coherence 
scale (see Enfilin 2003 for the related problem of the escape 
of cosmic rays out of radio cocoons). It is likely that, while 
AGNs and galaxies provide a potentially strong seed mag- 
netic field, there would still be a need for their subsequent 
amplification and maintenance against turbulent decay. 

Altogether, we adopt B BCC d = 10~ 7 G as a plausible es- 
timate of the seed magnetic field in the intracluster gas. 



APPENDIX B: THE CORRELATION 
FUNCTION OF THE FARADAY ROTATION 
MEASURE 

In order to calculate the autocorrelation function of the 
Faraday rotation measure, defined in Eq. 126H . we introduce 
coordinates {x,y, z) with the z-axis directed towards the ob- 
server, and those in the plane of the sky, X — (X, Y). We 
assume the magnetic field to be an isotropic, homogeneous, 
random field with zero mean value. Then its equal-time, two- 
point correlation tensor has the form (Bi(x,t)Bj(y,t)) = 
Mij(r,t), where 

M » = " rj ^f) Mtt{r,t) + ^M L (r,t) . 

Here (• • •) denotes the ensemble average, r — \x — y\, and 
Ti = Xi — yi (see Section 34 of Landau & Lifshitz 1975; Monin 
& Yaglom 1975). The functions M L (r,t) and M N (r,t) are 
known as the longitudinal and transverse correlation func- 
tion of the magnetic field, respectively. Since V • B = 0, 

We also assume for simplicity that the electron density is 
uncorrelated with the magnetic field and also constant over 
the field correlation length. This is consistent with the fact 
that random gas motions in galaxy clusters are quite sub- 
sonic. The correlation function of RM is then 

C(R) = (RM(Xi)RM(X 2 )) 

= K 2 n 2 c / / (B z {X 1 ,z 1 )B z (X 2 ,z 2 )}dz l dz 2 
Jo Jo 

= K 2 n 2 L J M zz {R,Qdt 
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Here we have assumed that L is much greater than the cor- 
relation length of the magnetic field, £ = Z\ — zi, R = 
\Xi-X 2 \ and r 2 = R 2 + ( 2 . 

For the sake of illustration, consider the longitudinal 
correlation function of the form 

Af L = i<B 2 >e X p 

which corresponds to the one-dimensional magnetic spec- 
trum of the form M k oc k 4 exp (-fc 2 Z|/2) (Monin & Ya- 
glom 1975). We note that Mk attains maximum at a 
wavenumber fc m = 2/is, (or a scale 2n/k m — ttIb), 
whereas the longitudinal correlation scale is given by ?l = 
t M L (0)] - 1 /" M L (r) dr = l B ^/2. 

Straightforward calculation then yields 

V% T s2 2 / r,2\ T i ( R_ 

\ 

where 

l-L/(V2l B ) 

(1 - s 2 ) exp (-s 2 / 2 ) « 0.88 , 

L/(V2l B ) 

with the numerical value obtained for L/(\^21b) 2> 1. 

The r.m.s. value of RM can be obtained from Eq. <BH 
or lHJl at R = 0: 

CRM = K 2 n c L J M N (R,z)\ R=0 d( 

= ^cK 2 n 2 (B 2 }Ll B , (B3) 
which is similar to Eqs 1271 and 1281 . 



C(R) = ^-cK"ni{B')Ll B exp I ] , (B2) 



